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Abstract. We study the Ginzburg-Landau model of superconductivity in three dimensions and 
for strong external magnetic fields. For magnetic field strengths above the phenomenologically 
defined second critical field it is known from Physics that superconductivity should be essentially 
restricted to a region near the boundary. We prove that the expected region does indeed carry 
superconductivity. Furthermore, we give precise energy estimates valid also in the regime around 
the second critical field which display the transition from bulk superconductivity to surface 
superconductivity. 



1. Introduction and main results 

The phenomenological Ginzburg-Landau theory of superconductivity successfully describes 
the behavior of a superconductor subject to an external magnetic field. Also, as pointed out in 
the celebrated work of Abrikosov, this theory predicted the existence of type II superconductors 
before they had been empirically realized, see |15| for a review of this physical topic for which 
A. Abrikosov was awarded the Nobel Prize. 

In the Ginzburg-Landau theory, the superconducting state of a sample is described by a 
complex-valued wave function ip and a vector field (magnetic potential) A such that the pair 
(■0, A) is a critical point of a specific energy (see (jl.ip below). The physical interpretation of 
ijj and A is explained by the microscopic Bardeen-Cooper-Schrieffer (BCS) theory as follows: 
is proportional to the density of superconducting particles and curl A measures the induced 
magnetic filed inside the sample. The rigorous mathematical justification of the connection 
between the Ginzburg-Landau and the BCS theory has only been established recently in |14j . 

The behavior of a type II superconductor is distinguished by three critical values the intensity 
of the applied magnetic field can have, that we denote by , and Hc-^ . These critical fields 
may be described in terms of the wave function ■0 as follows. Suppose H is the intensity of the 
external magnetic field applied to the sample. If H < Hqi , the material is in the superconducting 
phase, which corresponds to |0| > everywhere. If Hc\ < H < Hc2, the magnetic field 
penetrates the sample in quantized vortices (corresponding to zeros of ^). If < H < He,,, 
superconductivity is confined to (part of) the surface of the sample (corresponding to 1-01 very 
small in the bulk). Finally, if if > Hc^, superconductivity is lost, which is reflected by ■0 = 
everywhere. 

In the last two decades, much progress has been made in order to establish the aforementioned 
behavior of Type II superconductors by studying minimizers of the Ginzburg-Landau energy. The 
monograph |24| and references therein contains an analysis of vortices and the critical field Hc^^ ■ 
Concerning the analysis of the critical fields and He,, we mention |13| [8] (and references 
therein) . As one can see in [51 [23] , the Ginzburg-Landau model has a rich mathematical structure 
whose analysis requires a diversity of methods, and many of them have been developed especially 
for the study of the model. 

While a detailed study of the Ginzburg-Landau model in a two dimensional domain has been 
the subject of numerous papers, the study of the model in a three dimensional domain is much 
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less developed. Among the important problems that are still open in 3D is the calculation of 
the critical field Hc^ for general domains0. Also, when comparing with the existing results for 
2D domains, a precise localization of the wave- function if) is absent when the external applied 
magnetic field varies from Hc^ up to Hc^. However, for both 2D and 3D domains, a sharp 
characterization of the critical field Hc^ is given in |10| . In this paper, together with |12| I19| . 
we give a detailed description of the behavior of the wave-function tjj and its energy for external 
magnetic fields varying in the range above Hc^^ and up to Hc^ ■ The results in this paper concern 
surface (3D) superconductivity and the transition that happens close to from bulk to surface 
superconductivity. 

Let 17 C M'^ be a bounded and open set with smooth boundary which models a superconducting 
sample subject to an applied external magnetic field. The energy of the sample is given by the 
Ginzburg-Landau functional. 



2 

dx 



f3D(V,A)=£:3D(V.,A)= / [|(V-mi?A)Vp-AcVl' + 

Jo. L 

+ k^H'^ I \cmlA- dx . (1.1) 

Here k and H are two positive parameters whose physical interpretation are as follows, the 
number k is a material parameter, and the number H is the intensity of a constant magnetic field 
externally applied to the sample. As explained earlier, ip is a wave function (order parameter) 
and A is the induced magnetic potential. We take ip and A in convenient spaces as follows, 

where i!rjj^p(R'^) is the space introduced in (|1.2p below. Finally, f3 is the profile of the external 
magnetic field that we choose constant, (3 = (0,0, 1). 

Let i7^(M'^) be the homogeneous Sobolev space, i.e. the closure of C^(M^) under the norm 
^ '"^ ll^llij'i(lR3) II^^IIl2(k3). Let further F{x) = (— X2/2, xi/2, 0). Clearly divF = 0. 

We define the space, 

-^div,F(l^^) = {A : divA = 0, and A - F £ (R-^)} . (1.2) 
Critical points {tp, A) € H^{^}; C) x i^T^j^ p(]R^) of satisfy the Ginzburg-Landau equations, 
-(V -ZKHAftp = K^{l-\ip\^)ij in n, 

curl2A = -iKHA)ip)ln in M^^ (1.3) 

I N -{V - inHA)^ = on dO. , 

where 1^ is the characteristic function of the domain and N is the interior unit normal vector 
of dVt. 

For a solution (-0, A) of (jl.3p . the function ip describes the superconducting properties of the 
material and H curl A gives the induced magnetic field. 

The important class materials called Type II superconductors corresponds mathematically to 
the limit k — )• 00, see [HI I24| . 

We define the ground state energy, 

Eg.st(K,i7)= inf ^:^°(^,A). (1.4) 

(^,A)GHl(n;C)xHii^ j,(IR3) 

The leading order asymptotics of the ground state energy involves a function g constructed 
in |12) (the definition will be recalled in Section [2]) satisfying that g : [0, 00) — t- [—1/2,0] is 



-'^In a recent paper |18) . compactness results for the 3D functional valid in general domain are obtained, which 
allow to obtain the leading order term of Hci ■ Earlier results include a candidate for the expression of Hci in 
the case of the ball [5], and an expression of Hc^ in 'thin' shell domains in [B]. 
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continuous, increasing and there is a constant E2 < with 

5(6) = 0, V6>1, (1.5) 

g{b) = E2{1 - bf[l + o{l)], as 6 /I. (1.6) 

From |12) we have the following general (bulk) result as long as H/k is bounded from below 

\Eg.st{K,H)-gi§)\n\^\<CK^/\ (1.7) 

However, when H/k > 1, g(H/K) = and (jl.7p only gives a somewhat weak estimate of the 
energy. The value H = k corresponds to the phenomenologically described critical field 
mentioned previously. So one expects that around this value superconductivity should become 
concentrated near the boundary of the domain Vt and this should be reflected in the energy 
asymptotics. 

In this paper we will complete the study of the energy asymptotics for high magnetic field 
initiated in |12| by 

• giving the leading order energy asymptotics for H > k, which will be a surface energy of 
order of magnitude k, 

• studying the transition from 'bulk' to 'surface' dominated energy at i7 ~ k. We will see 
that this transition takes place at — k of order ^/k in the sense that 

— li H < K — /(k), where f{K)/^/K — )■ +00 as k — )■ 00, then the leading contribution 
to the energy comes from the bulk. 

— li H > K — /(k), where \imsupf._^^ /{k)/ ^/k < then the leading contribution to 
the energy comes from the surface. 

— li H = K — a^/H for some constant o > 0, then the bulk and surface contributions 
of the energy have the same order of magnitude. 

We now state the main results of the paper, Theorems 11.11 and 11.21 These require the intro- 
duction of some notation. 

• If X is a point on the boundary of Jl, then j^(x) denotes the angle in [0, 7r/2] between the 
vector /3 = (0, 0, 1) and the tangent plane to dVl at the point x. 

• If E [0,7r/2], C,{y') is the lowest eigenvalue of a magnetic Schrodinger operator in the 
half-space, see Section [3^ The function ^ is a continuous and strictly increasing bijection 
from [0, 7r/2] to [Gq, 1] where ©0 ~ 0.59 is a universal constant (the definition of Gq will 
be recalled in p.ip below). 

• If 1/ G [0,7r/2] and b < 1, the constant £'(b,z^) will be introduced in Section 13.3.31 
£^(b,i/) depends continuously on both b and vanishes when (^(z^) > b, and £'(b,i^) < 
otherwise. 

Below is a statement of the main result concerning the asymptotic behavior of the ground 
state energy. 

Theorem 1.1. Let [i : — ^ M 6e a function satisfying lim /x(k) = 0. Then there exists a 
positive constant kq, and a function err : — > M such that lim err(K) = and the following is 
true. If K > kq and H > k — ^{k)k , then the ground state energy in (jl.4p satisfies, 

Eg.si(K, H) = VkH / E (b, i^{x)) da{x) + E2\9.\ [k - H]\ + err(K) max (k, [k - H]\) . 
Jdn 

Here b = min [k/H, l), and da{x) is the surface measure on the boundary ofQ. 

The next theorem concerns the behavior of order parameters. We use the convention that a 
subset D C is smooth if there exists an open subset D C M.^ having a smooth boundary and 
such that D = Dnn. 
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Theorem 1.2. Let fx : M+ — )• M 6e a function satisfying lim /u(k) 
Then there exist a positive constant kq, and a function err : M_|_ - 
and if K> Kq and H > k — ^{k)k , then the following is true. 

(1) Ifii',A) G H^{n;C) X Hl^^piM.^) is a solution of ([LSI), then, 



and let D C 0, be smooth. 

1 such that lim err(K) = 0, 



1 



1^ dx < —K 



D 



^ ' _^(l7'^(^)) da{x)-E2\D\ 
1^ JdnnD 



K 

H 



1 



(2) Ifi^l^,A)eHHn■,i 



X HLf 



+ errf/t) max ( — 
is a minimizer of (jl.ip . then, 



' dx 



D 



annD 



E2\D\ 



K 

H 



+ err(«;) max ( — 



K 



[1.9) 



Remark 1.3. [Triviality for very large magnetic fields] 

Many previous works have addressed the question of the third critical field, see [8] and the 
references therein for a detailed review of this topic. We define 



(1.10) 



Hc.^{k) = inf {iJ > 0, for all H' > H the GL equations (fOD 
have only trivial solutions}. 

Here a solution is trivial if ■0 = 0. It follows from [8] that 

Hc^{k) = ©0 + o{k), as k — )• oo. 

Here, as we mentioned earlier, ©q ~ 0.59 is a universal constant whose definition will be recalled 
in ()3.ip below. So we can conclude that the ratio H/ k when k > 1 is always uniformly bounded 
from above in the regime where non-trivial solutions of (|1.3p exist. 

It was realized early in the Physics literature (see for example |15| Chapter 6.6]) that when 
studying the 3-dimensional Ginzburg-Landau model between and Hc^ there is a geometric 
boundary phenomenon that does not appear in 2D. When b = k/H ^ (©q, 1) (and k sufficiently 
large) superconductivity will in the 2D case exist on the entire boundary. But in the 3D model 
only a certain part of the boundary will carry superconductivity. This part can be described 
using the angle I'^x) (between the magnetic field and the tangent plane at the boundary point x) 
and the spectral function A local and linearized calculation suggests that superconductivity 
should only be present near the boundary section 

r{b) := {x e dn : C{i^{x)) <b}. (1.11) 

In previous works, notably |20j . the technique of Agmon estimates was used to prove that the 
superconducting order parameter ■0 will indeed decay rapidly away from T{b). Our results give 
the opposite direction, namely exhibit through Theorem 1 1 . 2 1 that minimizers are indeed nonzero 
(in an L^-sense) exactly in the vicinity of r(b). In order to realize this, notice that [k/ H —1]^ = 
by the condition on b and that E(b, z^) = if and only if b < C('^)- 

2. The universal constant E2 

Given a constant b > and an open set P C K^, we define the following Ginzburg-Landau 
energy, 



Gviu) = J (6|(V-fAo)up-|n|2 + i|u|^) dx . 



(2.1) 
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Here Aq is the canonical magnetic potential, 

Ao(2;i,X2) = ^(-X2,a;i) , V x = (xi, X2) G . (2.2) 

Given i? > 0, we denote by = {—R/2, R/2) x {—R/2, R/2) a square of side length R. Let, 

mo(fe,i?)= inf GkM. (2.3) 

u(iHl{KR;C) 

It is proved in |12| that 

• For all > 0, there exists a constant g{h) such that, 

i™ M = 9ib) ■ 

R^oo R'^ 

• The function , has a limit as 6 — > 1_, 
and -^<E2< 0. 

3. Reduced Ginzburg-Landau energy 

3.1. Harmonic oscillator on the half-axis. We denote by M+ = {i G M : t > 0}. For each 
real number ^ G M, we consider the harmonic oscillator 

H{0 = -^ + it-Cf in L\R+), 

with Neumann boundary condition, u'{0) = 0. Let ^i(^) denotes the first eigenvalue of H{£^). 
We define the universal constant, 

Go = inf^i(e). (3.1) 



The constant Qq satisfies the following properties [71 117): 

^ < 60 < 1, and 60 = Ail (Co), where = \/©o • (3.2) 
We introduce the function (fo G L^(M^) as follows, 

/ \ipoit)\'' dt = 1 , -<^'o'(t) + (t - eo)Vo(i) = 0o^o(t) inM+, ^o{0) = . (3.3) 

3.2. Magnetic Schrodinger operator in W^. We denote by = {{xi,X2,X3) G : xi > 
0}. Let v G [0, 7r/2]. Consider the magnetic potential 

E,y(x) = (0, 0, cos xi + sin 1/ X2) , x = (xi, X2, X3) G . (3-4) 

Notice that curl Ejy = (sin i/, — cos i^, 0) is constant. Geometrically, 1/ measures the angle between 
curl El, and the boundary of Mi]_. Consider the Schrodinger operator with constant magnetic 
field, 

C{i^) = -(V -iE^f in l2(]r3), (3.5) 

with domain 

D{C{v)) = {li G L2(IR3 ) . (V _ iY^^yu G L2(M+) , j = 1, 2 , = on {0} x M^} , 

We denote by Ci.^) the bottom of the spectrum of C{v): 

C{v) = m{a{C{v)) . (3.6) 
We collect below some properties concerning C(i/) (see e.g. [SJ Lemmas 7.2.1 & 7.2.2]). 

Lemma 3.1. Let Qq he the universal constant introduced in (|3.ip above. 

(1) The function [0,7r/2] 3 u ^ C{^) is monotone increasing, ("(0) = Gq and ^{tt/2) = 1. 
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(2) For all v G [0,7r/2), < 1. 

We denote by = (0, oo) x M. In connection with the analysis of the operator C{y), we 
introduce the two dimensional operator 

L{v) = -dl^-dl^ + {cosvxi + smvx2f in L^(R\), (3.7) 

whose domain is 

D{L{v)) = {n G F^(M+) : (cosi/xi +sinz^X2)-'n G L^{"^\), j = 1,2, d^^u = on {0} x M} . 
The link between the spectra of C{i') and Liu) is given below [8]. 

Lemma 3.2. Suppose v G (0, 7r/2). T/ien, 

(1) a(/:(z.))=a(L(^.)); 

(2) aess{L{v)) = [l,oo). 

Remark 3.3. Suppose that z/ G (0,7r/2). It results from Lemma [3. 21 and the properties of C(z^) in 
Lemma 13.11 that 

(1) Q{u) is an eigenvalue of finite multiplicity of L{v) ; 

(2) C(i/) is the lowest eigenvalue of L{v). 

Consequently, we can select a real- valued eigenfunction (f) G L^(R^) such that: 

(/>>0, / |(/>p(irE=l, / (|V0p + Kcosz^xi +sini/j;2)(/>n da:: = C(i^)- 

Moreover, using the technique of 'Agmon estimates' (cf. |17|). it is proved in |5l Theorem 1.1] 
that the eigenfunction (p decays exponentially at infinity as follows. If z/ G (0, 7r/2) and a G 
(0, y/l — C(^))) then there exists a positive constant C^^a such that 

J ^ e"\/^?+^(^|a^,0|2 + |a^.,</,|2 + Kcosi/xi + sini/X2)(/'|^ + dxidx2 < C^,a- (3.8) 

3.3. Reduced Ginzburg-Landau energy. Let u G [0,7r/2] and be the magnetic potential 
introduced in p.4p . For each £ > 0, we introduce the domains, 

K, = {-£,£) x{-i,i), U, = {0,oo)xK,, (3.9) 

and the space, 

Si = {ueL'^iUi) : {V -iB^)u£L'^{Ui), u = on {0,oo)xdKi}. (3.10) 
Let b G (0, 1] be a given constant. If u G Si, we define the Ginzburg-Landau functional, 

GbM'^) = (|(V - iE^)n|2 - b|u|2 + ^|n|4^ dx . (3.11) 

Associated with Qis,i,;t is the ground state energy, 

d{b,v-l) = inf gb,^;Kn). (3.12) 

3.3.1. Preliminary properties of minimizers. 

Lemma 3.4. For eac/i G [0, 7r/2], /ei C(z^) 6e as defined in ()3.6p . Lei £ > be given. 
Ifi>< C(i^), i/ien 

(i(b,i/;£) = 0, 

where d{b,i';i) is the ground state energy introduced in (j3.12p . 

Proof. By using ti = as a test function, it is clear that d(b,i/;^) < 0. On the other hand, the 
min-max principle and the condition on b show that d{b, v^t) > 0. □ 

Remark 3.5. It results from Lemma[3]4]and the properties of ^(i^) in Lemma[3TT]that d{b, v]tj =0 
in each of the following cases: 
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(1) b < Go and z> G [0,vr/2] ; 

(2) b £ [Go,!] and z^ = 7r/2. 

Theorem 3.6. Let Go be the constant introduced in (|3.ip . Suppose that b G [Go, 1] «s a given 
constant. 

For all V G [0, 7r/2] and £ > 1, the functional Gb^ul in (|3.1ip admits a minimizer fb^^-£ £ Se 
satisfying, 

Gb,y/{^b,v;l) = d{b, v; t) , \\ipb^y-i\\L°-=-(Ue) < 1 • (3.13) 
Furthermore, there exists a universal constant C > such that, if v £ [0, 7r/2] and I > Q, the 
minimizer (pb,u;e satisfies, 



I 

JU, 



Xi 



luen{xi>A} (lnxi)2 
Proof. Let m > and consider the energy functional, 



(I (V - iE^)ipb,u;£\^ + \ipb,u;i\'^ + xl\ipb,u-l\^) dx < Cf . (3.14) 



{u)= f ( |(V-iE^)np-b|up + -|un fix, (3.15) 



where is the magnetic potential in (j3.4p and, 

Ui^m = (0, m) X , = (-^, £) X {-£, I) . 
We define the ground state energy, 

t{i,m) = inf Ge,miu) , (3.16) 

where 

Si^„i = {n G H^iUi^m) : n = on (0, m) x and {?n} x iT^} . (3.17) 

The proof of Theorem 13.61 consists of showing that t{i, m) — )• d{b, v]li) and that a minimizer (yJ^^m 
of Qt^rn converges to a minimizer of Qit,u;t as m — )■ oo. That will be done in several steps. 

Step 1. In this step, we prove that t{i,m) — )• d{b,v;i) as m — )• oo, where d{b,i';£) is the 
ground state energy introduced in (|3.12p . Recall the space Sg introduced in ()3.10p . Since every 
u G Se^m can be extended by to a function u G Sg, we get that, t{i, m) > d{b, v; £). So, we need 
only prove that, 

limsupt(^,m) < d(b,z^;£) . (3.18) 

m—¥oo 

Let (ipn) C 5^ be a minimizing sequence of Gb,u;£: i-^- 

d{b,i';£)= lim Gb,u;£{Vn) , 

n— >oo 

where Gb,u;e{fn) is the functional introduced in (|3.1ip . 
Let r] G C^(M) be a cut-off function such that, 

0<7?<linIR, ?7 = 1 in [-1/2, 1/2] , supp C [-1, 1] . 

Let Tjmixi) = r/(xi/m), then r]m{xi)ipn{x) G Si^m and consequently, 

t{£,m) < Gb,v;l{'^m'^n) ■ (3.19) 

Using the inequality, 

|(V - iE^)(7?^(^„)|2 < (1 + e)|7?„(V - iE^)¥^„|2 + 2e"V^ml^l¥'n|^ 
valid for all e G (0, 1), together with the fact that < r]m < 1, we get the following estimate, 

t{i,m)<{l + e)Gb,,.A'Pn) + ^^\W\\lo.m[ \ipn\'^dx + b [ {1 - r]i + e) \ipn\'^ dx . (3.20) 
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Taking limsup on both sides of (|3.20p . we get (by using in particular Lebesgue's dominated 
convergence theorem), 

limsup t{£, m) < {1 + e)gh,uA^n) + be / \(pn\'^ dx . 

m— >-oo Ju^ 

Taking successively e — )• 0+ then ?i — )• oo, we get the estimate in (|3.18p . 

Step 2. Since Ui^m. is a bounded domain, it is obvious that Qt^m has a minimizer i-pi^m^ see 
e.g. [SI Chapter 11]. The function (^^^m satisfies the following equation, 

- (V - ?E,,)^V3£^m = b(l - , inf/^^m. (3.21) 

A simple application of the maximum principle yields ^ 1 everywhere. We will prove the 

following estimate, 

I TT^Nl (K^ - ^^y)^t.m? + + dx < Ce , (3.22) 

^c/£,,nn{xi>4} (ln2;i)^ 

valid for all £ > 0, where C is a universal constant. 

We obtain the estimate (|3.22p by following a construction similar to |21| . Consider a function 
X G C°°(M) such that suppx C (0, oo). Using the equation in (|3.2ip and an integration by parts, 
we get. 



/ {\{"^ - i'^u)XVt,rn? -^XVt,rn? + bX^Wt,rn\'^) dx = 1 W?Wt,m? dx . (3.23) 

J Up yyi J Up ^ 

(K 

magnetic field in L^(M'^) is equal to 1, we get the lower bound. 



Since xVlrn ^ H^(R^) and the first eigenvalue of the Schrodinger operator with constant unit 



/ \{V - iB^)x^i,m\^ dx > / 



\X^i^rn? dx. 



Inserting this lower bound into (|3.23p and remembering that b G [©o, 1], we deduce the following 
estimate, 

X^\f£,m\'^dx<^f \x'\^\<fi,m\^ dx . (3.24) 
We select the function x so that, 

3/2 

X(xi) = if xi < 0, x{xi) = P—i^xi>A. (3.25) 

mxi 

The function x consequently satisfies 

3 /x7 

< x'ixi) < for ah xi > 4 . (3.26) 

21n2;i 

Using the properties in (|3.25p and (j3.26p of the function Xi we infer from (|3.24p . 

^3 

1 I |4 J 

7KWe,7n\ dx 



1 

< - 
- b 



9 /■ XI . .2 , f 

7 / n '^m,rn\ dx+ / 

4 Jc/,,„n{xi>4} (inxij Jue,„,n{xi 



<4} 

1/2 



\x'f\^£,m\'^dx 



TT^IVf.mr^^ e + Ci', (3.27) 
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where C is a universal constant. We explain how we get the estimate in ()3.27p . Actually, using 
that \'ft,m\ < 1, we get, 

/ W?Wi,m?dx < ||x'IIl-'([o,4])(16^^) • 

Jf/f,„n{a;i<4} 

Also, using a Cauchy-Schwarz inequality, we have, 
/ 7r—^-T^\^e,m\'^dx 

<([ 1 A'^'ff 4 i4,,y'' 

~ \Jut,^n{xi>i} xiilnxiY J \y[/,,^n{xi>4} (lnxi)2 J 

thereby yielding (|3.27p . 

As a consequence of (|3.27p . we deduce the following estimate, 

V£>1, / {\ipi,m\^ + xi\ipe,m\^) dx < C£\ (3.28) 

for a possibly new universal constant C. Thus, to finish the proof of p.22p . we only need to 
prove that 

/ TT^^li"^ - i'Eu)^i,m.\^ dx < Ce\ (3.29) 

JUi,mn{xi>4:} (hixi)^ 

To get (|3.29p . we select x so that ^(^^i) = y/^/^^xi if xi > 4, and ^(^^i) = if xi < 1. Using 
the estimate, 

/ x\V - dx < 2 [ {\{V -iB^)xiPi,m\'^ + \x'\^\v^i.m\^) dx, 

together with the estimates in (j3.23p . (|3.28p and the properties of we get easily the estimate 
in (13:29]) . 

Step 3. In this step, we prove that there exists a minimizer ipi of Qi,,u\i satisfying the estimate 
in (imH . 

Starting from the bound |(/J^,m| ^ 1 and the equation ()3.2ip . we get by a compactness and a 
diagonal sequence argument that there exists a function ^pi E C^{Ui) such that, as m — )■ oo, 

^Lm^^i in C^(K), 
for any compact set K dUg. Furthermore, ipg satisfies the equation, 

-(V-iE^)V£ = Kl-|'/5£|^)'/5^, int/^, (3.30) 
together with the boundary conditions, 

^— ^ = on {0} X K'^ , = on (0, oo) x dK^ . 

Using the estimate in p.22p together with the bound |(/7£,mJ < we get. 



/ (|(V - iB^)^i,m? + \^l,m?) dx < Ce 
•J Up ^ 



Passing to a subsequence, we conclude that, as m — )■ oo, 93^^^ ^ '■Pt ^ind (V — lE^)^?^^, 
(V — iE^)ip£ weakly in Lp'iUi). Consequently, we get that (yj^ is in the space 5^, 



Wl,m\^ dx ^ I dx as m — )■ 00 , 



and (pi satisfies the estimate in (|3.14p . 
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Using the equation (|3.2ip and an integration by parts, it is easy to check that, 

t{e,m) = I \fe^m\'^dx. 
Letting m — t- oo and using Step 1 above, we get, 

d{b,i^;i) = -7; [ I'^el'^dx. 
But, thanks to ()3.30p . an integration by parts yields, 

thereby obtaining that (pi is a minimizer of Qb^u/- D 
3.3.2. A rough estimate of the ground state energy. 

Lemma 3.7. There exists a universal constant Ci > such that, for all b £ [Qq, 1], u € [0,7r/2] 
and i > 0, we have, 

d{b, v; e) > Ci min(C(i^) - b, O)^^ . 
Furthermore, if v ^ [0, 7r/2) satisfies C,{u) < b, there exist positive constants C2 and Iq such that, 
if i > £q, then, 

d{b,u;i) < -C2f 
Remark 3.8. It resuhs from Lemma [3.71 that 

-00 < mf -5 < sup -5 < , 

provided that b and 1/ satisfy the conditions in Lemma 13.71 
Proof of Lemma \3. 71 

Lower bound: Let ut, G 5^ be a minimizer of Qb^v;l- Since Uf, = on (0, 00) x dK^, we can 
extend u\, to a function u^, by letting 5[, = in (0,oo) x (M^ \ Kf). The function 5(, is in the 
form domain of the operator C{i') introduced in (|3.5p and, 

d{b,v;(.) = Qb,uA'^b) 

= j ^ ^|V - iE^)nbP - b|SbP + ^l^bl"^^ dx . 
Using the variational min-max principle, we get. 



d{b,v;l)>{Q{v)-b) \ub\^ dx = {({ly) - b) \ub\' dx . 

By Theorems we have / < where C, is a universal consU,,.^ If CM - i, < 0. 

then 

d{b,u;£) > Ci{C{i^)-b)f. 

If C(i^) — b > 0, then Lemma 13.41 tells us that d(b,i^;£) = 0, thereby obtaining the lower bound 
announced in Lemma 13.71 

Upper bound: We construct a test-configuration as follows. Suppose v G (0, 7r/2). Let 
X G C^i^) be a cut-off function such that 0<x<l, X = lin [-1/2,1/2] and x = in 
M\ (-1,1). For ah v = (s,t) G M^, we set 

miv) = Xiis)xeit), Xdx) = x{^) (xGR). (3.31) 

Let (p be the eigenfunction introduced in Remark 13. 3| and let M > be a given real number. 
For j G Z we set Cj = Mj. The number of such Cj that belong to the interval {—£,£) is given 
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by N{i,AI) = 2{i/M} + 1, where {r} denotes the largest integer strictly less than r G M. For 
X = (xi,X2,X3) E M^, we define our trial state fi{x) to be 



{e/M} 

ili{x2,X3) ^ exp(-zcj sin(i/)3;3)(^(xi,X2 - Cj 
j=-{e/M} 



We compute Gb,u;eitfe), where t > is a constant whose choice will be specified later. Actually, 



Gb,uA*M = ''^^ 'ne{x2,X3) ^ ex.p{-icj sm{u)x3) {dx^^){xi, X2 - Cj) 

j=-WM} 

WM} 

+ r]i{x2,xs) ^ exp{-icj sm{u)x3) {dx2(p)ixi, X2 - Cj) 
j=-{e/M} 

{l/M} 

+ dx2{ili){x2,X3) ^ exp[-icjsm{u)x3)(j){xi,X2 - 
j=-{e/M} 

WM} 



+ 



m{x2,X3) ^ exp(-zcj sin(i/)x3) (cos(z^)2;i + sm.{v){x2 - Cj))(j){xi,X2 - cj] 
j=-WM} 

WM} 

+ dx^[r]g{x2,X3)) ^ exp[-icjsm{u)x3)<j){xi,X2 - Cj] 

j=-{i/M} 
WM} 

'niix2,X3) ^ exp(-icj sin(i/)x3)(^(xi,X2 - Cj 

j=-WM} 

H — t^%(x2,X3) exp(— zcj sin(i/)x3)(;i(xi, X2 — Cj) dx . (3.32) 

j=-{i/M} 



We define Graainiifi) and Grestitfi) as foUows, 



WM} 



Qma.initfi) = t"^ \rilix2,X3)f ^ |^| (c?a:^ (/)) (xi , X2 - Cj) ^ + | (9^2 (xi , X2 - 

j=-{i/M} 



+ |(cos(i/)xi +sin(i/)(x2 - Cj))(/)(xi,X2 - Cj)\ 



b\(j){xi,X2 - Cj)P + -t^|0(xi,X2 - Cj)\'^ 



dx 



and 



Grest{tfe) = Gb,u;e{tfe) — Gmain{tfe)- 



We will show that Qrestitfe) is small compared to Qmainitfe). We start by obtaining an upper 
bound on ^main(i/^)- Using that %(x2,X3) = Xe{x2)xe{x3) and < Xe{x2) < 1 we get for a fixed 
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integer j 



|%(2;2, 2:3)1^ \{dxi(p){xi,X2- Cj)f + \{dx2(l))ixi,X2- Cj)[ 



+ \ [cOs{v)xi +Sm{v){x2 - Cj))4){xi,X2 - Cj)|^ + -t^\4){xi,X2 - Cj 
<t^ \Xiix3)\'^ dX3 \idxi(p){xi,X2 - Cj)f + \{dx2(l))ixi,X2 - Cj)\ 



dx 

2 



+ \ [cos{u)xi + sin(z^)(rc2 - Cj))(/){xi,X2 - Cj)^ + X2 - Cj)]"^ dxidxQ 



(xi, 3:2)1 dxidx2 



where A = / dz. For the negative term, we have 

Jr 



-bt'' I \'ne{x2,X3)(j){xi,X2 - Cj)\ dx 

lUi 



-bt' 



-bt' 



\r]£{x2,X3)(l){xi,X2 - Cj)\'^dx 



Xeix3)\'^dxs / \x£ix2)(l){xi,X2 - Cj)\'^ dxidx2 



= -b£Xr + biXr / {1 - \xe{x2)\ )\(t>{xi,X2 - Cj)\^ dxidx2 

Using the exponential decay of (p we find that for ah s > 

-bf I \'ni{x2,X3)4>{xi,X2-Cj)\^dx = -beXt"^ + o{l/i') 



as ^ —7- 00. We conclude that 

Gm.in{tfe) < N{£, M)£Xt^ [c(i^) " ^' + 
as £ — ?■ 00. We select t sufficiently small so that 



{xi,X2)\ dxidx2 



+ N{£,M)o{l/e) 



C{u) - b + 



|(/>| dxidx2 < 



C{u) - b 



to get 



Gm^initfe) < ^Nii, M)£Xt\C{u) -b)+ Nii, M)o{l/f) < ^eXt\Q{v) - b) (3.33) 
as £ — )• 00. 

All terms in Grest{tfe) can be taken care of in the same way, using (|3.8p (except for the term 
with dx2'neix2,X3) which is exponentially small due to the decay of (p). For this reason we only 
show how to handle the terms in Grestitfe) that comes from the first absolute value in (|3.32p . We 
meet terms like 



I |r/Kx2,X3)pe-^<^^-^'^('^)^3(a^^^)(^^^^2_c.)e-'=-W3(9^^ (3.34) 

where j ^ k are such that cj and Ck both belong to {—i, i). We assume that j < k and estimate 
the absolute value of this integral by 



t Xi \{dxi)(l){xi,X2 - Cj){dxi)(l){xi,X2 - Ck)\dxidx 



2- 
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Next we decompose the integral into two terms, one where X2 < {cj + Cfc)/2 and another one 
where X2 > {cj + Cfc)/2. Using p.Sp we get 

j(xuX2)ml l(^a;i)<^(a:i,a;2 - Cj){d^^)(j){xi,X2 -Ck)\dxidx2 

X2<(cj+Ck)/2 



Here we have used the L°° bound on the first exponential, and then the Cauchy-Schwarz in- 
equality together with the exponential decay on from (|3.8p . The same bound is true for the 
integral where X2 > {cj + Cfc)/2. Looking at the case j > k we get the same bound, but with 
j — k instead of fe — j in the exponential. We thus find that ()3.34p in absolute value is bounded 
by 

We want to sum over all integers j ^ k running from —{£/M} to {£/M}. For a real number 
a < 1, the sum 

2 Bn - j)a^ = (a" - 1 + n(l - a)) < 

Our sum transforms into this sum with a = e~"^^/^ and n = {i/M} < IjM we find that the 
error coming from the first absolute value in (|3.32p is bounded in absolute value by 

-qA//4 



' M{l-e-»^I^Y' 



Doing the same estimates for the other terms, we find that there exists a /3 > and a constant 
C > (not depending on M) such that 

Qrest{tft)<Ct'\e-P^e. (3.35) 

Combining (|3.33p and ()3.35p we find that 



Choosing M large enough we find that there exist Iq and C2 > such that ii i > £0 then 
Qb,u;e{tfe) ^ —6*2^^. This finishes the proof for u £ (0, 7r/2). 

Suppose that = 0. In this specific case, we select the trial function as follows, 

fi{x) = m{x2,X3)ipo{xi) , 

where (po is the ground state introduced in (j3.3p , and rji is the cut-off function introduced in (j3.3ip . 
Consider t > 0. An easy computation yields, 

for some constant C. Selecting t sufficiently small so that Gq — bH / 1 990(3^1)!^ dxi < , 

2 Jr_^_ 2 
we get the upper bound that we want to prove. □ 
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3.3.3. The thermodynamic limit. The aim of this section is to prove the following theorem. 
Theorem 3.9. 

(1) Suppose b G [©Oil] o-'iT'd v G [0, 7r/2]. There exists a constant E{b,i') such that 

Furthermore, E{b,v) < if b > (^{ly) and E{b,i') = otherwise. 

(2) There exist positive universal constants Iq and C such that, 

E{b,iy)<^^^^^<E{b,u) + -^, Vi/G [0,7r/2], y £ > io . (3.36) 

The proof of Theorem 13.91 relies on the following abstract lemma, which is proved in |12j . 

Lemma 3.10. Consider a decreasing function d : (0,oo) — )• (— C)O,0] such that the function 
f : {0,oo) B e ^ ^ £R is bounded. 

Suppose that there exist constants C > and > such that the estimate 

f{ni) > /((I + a)e) -c(^a + -1^ ) , (3.37) 

holds true for all a G (0, 1), n G N and I > £q. 

Then f(£) has a limit A as £ —)• oo. Furthermore, for all i > 2£q, the following estimate holds 
true, 

2C 

+ (3-38) 

In order to use the result of Lemma [3. lUI for the energy d{b,v,£), we establish the estimate in 
Lemma 13.111 below. 

Lemma 3.11. Let b G (Go, 1) and v G (0, 7r/2) such that b > Ci'^)- There exist universal 
constants C > and £o ^ ^ such that, for all £ > £q, n G N and a G (0, 1), we have, 

d{b,u;ni) ^ d{b,iy;{l + a)£) C 



Proof. Let n > 2 be a natural number. If a G (0, 1) and j = (ji, j2) G we denote by 

Ka,j = Iji ^ I 32 ' 

where 

V m G Z , Im = ^"T- — n — 1 — a — (1 + a), m — n — 1 — a+(l + a) 
Consider a partition of unity {xj) of such that: 

(J 

^\Xj? = l, 0<Xi<l inM^, suppxi C -FCaj, |Vxil<-, 

3 

where C is a universal constant. We define X£j(^) = Xj{^/^)- Then we obtain a new partition 
of unity Xi,'j such that suppx£,j C /C^j, with 

/C^j = {£x : X G Ka,j] . 

Let J = {j = {ji-,32) G : 1 < ji,j2 < n} and = {—n£,n£) x {—n£,n£). Then the family 
i^e,j)jej is a covering of K^e, and is formed exactly of squares. 

We restrict the partition of unity {xi,j) to the set Knt = {—n£,n£) x {—n£,n£). Let Uni be a 
minimizer of ()3.1ip . i.e. Q\i^u]ni{uni) = d{b,v]n£). We have the following decomposition formula, 

Qb,u;nt{Uni) > ^ {Qb,v;ni{X£,jUni) " || \^ X^,3\Uni\\\2^^K„l)) ■ (3.39) 
3&J 
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By magnetic translation invariance, we get Qb,u;ni{Xi,jUnl) ^ d{b,v;{l + a)l). Therefore, it 
results from (|3.39p . 

d(b, v; nl) > n^d{b, v; (1 + a)i) - \\Vxi,j ^mWh^K^^) • (3.40) 

j 

Using Theorem 13.61 we get further, 

2 Cn^ 



d(b, V- ni) > n''d(b, v- (1 + a)l) 5- . (3.41) 



Therefore, 



d{b,v\nl) ^ d{b,v-{l + a)l) C 



□ 



Remark 3.12. Using the lower bound in Lemma 13.71 we infer from Lemma 13.111 the following 
lower bound, 

d{b,v-nl) ^ d{b,v,{l + a)^) 1 



[nlf - ((l + a)£)2 V a^^^ 

where the constant C > is universal. 

Proof of Theorem 



Let fU) = ^i^l^. Thanks to Lemmas 13.71 and 13.111 we know that the functions f{t) and 

d{t) := d{b, ly; £) satisfy the assumptions in Lemma [3 .101 Consequently, f(i) has a limit E(b, v; I) 
as ^ — )■ 00 and, 

for all I > Iq. Here C and £q are constants independent of b E [0o, 1] and S [0,-7r/2]. 

It remains to establish a lower bound of f{i)- Let n G N. By using a comparison argument 
and magnetic translation invariance, we have that 

d{b,iy;ni) < n'^d{b,iy;i) . 

Consequently, 

d{b,u;i) ^ d{b,u;ni) 
£2 - (n£)2 • 

Making n — t- 00 in both sides above, we conclude, 

'-^>E(>,,^). (3,42) 

Lemma 13.71 tells us that if b > C(^)) then d{b,i';i) < 0. In this case, we get i?(b, i/) < as 
consequence of (|3.42p . On the other hand. If b < C(^)) know by Lemma 13.41 that E{b,i') = 
0. □ 

3.3.4. Properties of the function E{b,i'). Theorem 13.91 provides us with a limiting constant 
E{b,i') G (—00,0] defined for b E [©0)1] and E [0, 7r/2]. In this section, we will study 
properties of E{b, v) as a function of b and v. 

Theorem 3.13. For all {b,v) E [©o,!] x [0,vr/2], let the constant E{b,v) he defined as in 
Theorem \3.9l Then E{b, u) satisfies the following properties. 

(1) Given v E [0, 7r/2], the function [©0, 1] 9 b 1— )■ Eifo^v) is continuous and monotone 
decreasing. 

(2) Given b E [©0, 1], the function [0, 7r/2] 9 1— )■ E{b,v) is continuous. 
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Proof. Continuity and monotonicity of b i— )■ E{b, v): 

Let b E [00, 1] and e E M such that b + e E [Qq, 1]. Recah the definition of the functional Gb^u/ 
in ()3.1ip together with the associated ground state energy d{b,i';£). It is easy to check that. 



dx , (3.43) 

'u, \ ^ J 

valid for any u E Sg. In particular, setting successively u = ^t,,u;i then u = in (|3.43|) and 
using the properties in Theorem 13.61 we get the following estimate, 

\d{b + e,u;i) -d{b,iy;£)\ < Cef , 

for some universal constant C. Remembering the definition of £'(•,£/), we get, 

\E{b + e,v)- E{b,v)\ < Ce, 

thereby proving the continuity of £'(•, v). 

To obtain monotonicity of E[b,v), we suppose that e < and we set u = ^Pb+e.u-/ in (|3.43p . 
Thanks to Theorem 13. 6| we know that |u| < 1 and consequently, 

d{b + e, u; i) > d{b, [ dx 

> d{b,u;£). 

Dividing both sides above by 4£^ then letting ^ — )• oo, we get, 

E{b + e,iy;£) > E(b,i/;£), 

which proves that E{-,v) is monotone decreasing. 
Continuity v ^ E{b,v): 

Let V E [0, 7r/2] and e E (—1,1) \ {0} such that + e E [0, 7r/2]. We want to prove that 
E{b, v + e) ^ E{b, u) as e 0. 

Recall the definition of the magnetic potential Ej, in ()3.4p . Notice that, if u E Si, we have the 
following estimate, 

\gb,u+sAu) - Ob,uAu)\ < \e\ [ \iV - iB^)u\'^ dx + 2\e~^\ [ |(E^+e - E^)^!^ . (3.44) 

JUi JUe 

Using the bounds, 

I cos(z^ + s) — COS i^l < e , I sin(z^ + e) — sin iy\ < e , 

we get, 

\E^+e{x) - 'E,Ax)\ <£{\xi\ + \X2\) , V X = {X1,X2,X3) E M.^. 

Consequently, we infer from ()3.44p the following bound, 

\Gb,u+eA'^) - Gb,uAu)\ < \£\ [ \{V - i'E^)u\'^ dx + 2\e\ [ (jxip + jxap) jnp dx . (3.45) 

JUi JUe 

Let C e ^^(M) be a cut-off function satisfying < C < 1 in M, suppC C [-1, 1] and C = 1 in 
[— ^, Let Q{x) = C{xi/£) and ipi,^i/+£;£ E 5^ be a minimizer of Gb,u+e-/ given by Theorem 13.61 
Applying p.45p with u = Qv^i ^^nd using the decay estimate of Theorem 13. 6| we get for all ^ > 1, 

Gb,.+eACm,u+eA > d{b, v; i) - , (3.46) 

where C > is a universal constant. Using the estimate in (j3.36p . we get for all £ > io, 

Qb,u+eA^m,u+eA > (4^')^(b, v) - C\e\£^ , (3.47) 
where £q is a universal constant. 

We estimate the term Qb,v+eAit-'-Pb,v+eA from above. Actually, an integration by parts give 
us the following identity. 



Qb,v+eA'^t-Vb,y+eA = [ \Ce\'^\^b,u+£;e\'^ dx + [ Cl(oC|-l 



\ipb,u+e;e\'^ dx . 
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Using the decay estimate of Theorem 13.61 we get, 

Gb.u+e-ACifb.u+e-/) < C - - [ Ifb-u+e'/l"^ dx . 
The Ginzburg-Landau equation satisfied by '^b.u+e\i gives us 

Gb,v+e;i{'^b,u+e-/) = --^ [ \V^b,u+e;e\'^ dx. 

^ JUt 

Consequently, we get the upper bound, 

Qb,u+e-/{Ct'Pb,u+e;d < d{h,V + e'l) + C . 
Using the estimate in p.36p . we get further, 

Qb,u+eACm,u+e;d < {^^^)E{b, U + e) + Ct''' . 
Inserting this upper bound into (|3.47p . we get, 

E(h, v^e)> E{b, u) - C\e\£^ - Cr^/^ . 
Taking successively liminf then lim on both sides above, we get, 

liminf £;(b,z/ + e) > E{b,v) . 

In a similar fashion, by applying (|3.44p with u = Cefb,u;e and Vb,iy;e being a minimizer of Qb,u;l, 
we can prove that limsup£^(b, v + e) < E{b, v). This gives us that 

\im E{b,i^ + e) = E{b,i^) , 
and thereby proves the continuity of the function E{b,-). □ 

4. Auxiliary results 

We collect some results that are necessary for controlling the errors resulting from various 
approximations. The first result is an estimate obtained in jJl Lemma 3.2 and Theorem 3.3] (for 
a different method see [HI Theorem 12.3.1]). 



Lemma 4.1. There exists a constant Ci > such that if (^, A) G H {Q) x i7^j^p(M") is a 
solution to (11.31). then 



Here, 

1 



A = max 



K IH 



(4.1) 



The next lemma is taken from [51 Lemma 10.33], which, together with Lemma [4. H give a good 
estimate of || curl(A — F)||j;^2( 



Lemma 4.2. There exists a constant C > such that, if (ipjA) G H^{^l;C) x if^;yp(M'^) is a 
solution of (jl.3p . then, 

||curl(A-F)||^2(K3)<g||V|li4(^,), (4.2) 

for all K > and H > 0. 

The next result is Theorem 14.31 below which is proved in |12| . Similar estimates to those in 
Theorem 14.31 are also given in |20) . 
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Theorem 4.3. Suppose that < Amin < Amax- There exist constants kq > 1 and Ci > such 
that, if 

^ /^O ) Ajniii < — < Ajnax i 

and A) G i?i($7;C) x if]j^ p(]R2) is a solution of ([O]), then 

||(V - mFA)^llp(^) < Ci^/^||^||ioo(o) , (4.3) 

||A - F||iy2,6(n) < Ci (^11 curl(A - F)||i2(iR3) + -^||V'|U6(n)||V'llL-(n)^ , (4.4) 
||A-F||pi.i/2(f^) < Ci (l|curl(A-F)||i2(iR3) + -^||V'||L6(f^)||V'||i-(f^)^ , (4.5) 

II curl(A - F) 11^0.1/2(11) < Ci (^11 curl(A - F)||^2(k3) + -i=||^||^6(f,) || V||l^(q)) • (4.6) 

Combining the results in Lemma 14. H Lemma 14.21 and Theorem 14.31 we get the following 
corollary. 

Corollary 4.4. Suppose that < Amin < Amax- There exist constants kq > 1 and Ci > such 
that, if 

K 

1^ ^ 1^0 1 Amin ^ ~ !^ Amax i 
H 

and (V', A) G i?^($7;C) x if]j^p(]R3) is a solution of ([LS]), then, 

l|A - F||f;l,l/2(f^) < Cl . 

5. Lower bound 

In this section we will prove Theorem 15.11 below . whose statement requires some notation. Let 
D C r2 be a given open set such that there exists a subset D of M'^ having smooth boundary and 
D = D r\VL. For all a > 0, we assign to D the following subset of ^2, 

Da = {x(^n : dist(x, D)<a}. (5.1) 

We introduce the following functional, 

So{u,A;D) = j (^|(V-iA)n|2-K2|n|2 + :!^|u|^^ dx, (5.2) 

where u G H^{Q; C) and A G i^^J^^ p(M'^). If D = $7, we omit the dependence on the domain and 
write £q{iP, A) for £o{ip, A; il.). 

Theorem 5.1. Suppose that the magnetic field H is a function of k such that, 

H H 
1 < lim inf — < lim sup — < oo . 

Let K 3 M+ I— 7- a{K,) G M+ be a function satisfying lim o^k) = 0. Then, for any solution 
(V',A) G H'^{n;C) X if]i^p(]R3) o/ ([H]) ancf any function h G C^(r2) satisfying ||/i||Loo(f^) < 1 
and supp /i C -Da{K) ; the following asymptotic lower bound holds, 

£o{h->P,A)>V^ _ E(—,u{x)) da{x) + E2\D\[K-H]\ + o max (k, [k - H]^) 
JDnan \H J \ J 

as K — )• oo . (5.3) 

Here da{x) is the surface measure on the boundary of Vt, E2 < ^ is the universal constant 
introduced in (|2.4p , and £q is the functional introduced in ()5.2p . 



3D GINZBURG-LANDAU FUNCTIONAL 



19 



The proof of Theorem 15.11 is a bit lengthy and is divided into several subsections. We will 
introduce two parameters 

a = a{K)£ (0,1), 6 = 5{k) £ (0,1) , 

that will be used along the proof. Different conditions on these parameters will arise so as to 
control the error terms correctly. The choice of the parameters will be fixed at the end of the 
proof (in ()5.18p below) so that they satisfy the aforementioned conditions, and are negligible in 
the limit of large k. 

Throughout the section, {^,A) G H^{^1) x if^j^p(M^) will always denote a solution to (|1.3p . 

5.1. Splitting into bulk and surface terms. We introduce smooth real- valued functions xi 
and X2 such that Xi + xi = 1 ^) 



Xiix) 



1, if dist(x,aJ7) < 5/2, 
0, if dist(x,a!^) > 5, 



and |Vxj| < C/6 for j = 1,2 and some constant C independent of S. Using the IMS decompo- 
sition formula and the fact that /q(Xj(^)^ ~ Xj{x)^)\'^\'^ dx > (since < Xji^) ^ 1)? we find 
that 

2 r 

£o{H,A.)>£o{xiH,A)+£o{X2H,A-)-y^ \VXj?\H? dx. (5.4) 

To estimate the error we use that ||V'||oo ^ 1 and the fact that the measure of the support of 
Vxj is bounded by C5 for some constant C, so get 

2 r 

V / \Vxj?\H?dx<C5-^. (5.5) 

Next, we estimate separately the terms £Q{xihil),A) (surface energy) and £o{x2hil^ , A) (bulk 
energy) . 

5.2. The surface energy. The estimate of the surface energy requires two steps, a decomposi- 
tion of the energy via a partition of unity, then passing to local boundary coordinates that allow 
us to compare with the model case of the half-space that is studied in Section |3l 

5.2.1. Boundary coordinates. We introduce a system of coordinates valid near a point of the 
boundary. These coordinates are used in [TB] and then in |22) in order to estimate the ground 
state energy of a magnetic Schrodinger operator with large magnetic field (or with small semi- 
classical parameter). 

Consider a point xq E dVl. After performing a translation, we may assume that the Cartesian 
coordinates of xq are all 0, i.e. xq = 0. Let y be a neighborhood of xq such that there exists 
local boundary coordinates {y2^yz) in VF = F n50, i.e. there exists an open subset U of M2 ^nd 
a diffeomorphism (f) : W ^ U, (p{x) = {y2, y^)- Denote by N the inward pointing normal at the 
point 4>~^{y2^yz) £ d^- We define the coordinate transformation <I> as 

{xi,x2,x-i) = ^~'^{yi,y2,yz) = <A~"^(y2,y3) + yiN. 
The standard Euclidean metric = dxj (8> dxj transforms to 



ao= 2^ gjkdyj ® dyk 



dys ^dy3+ ^ [Gjk{y2, ya) - 2yiKjfc(y2, ya) + yi^i/c(y2, ys) 



2<i,fc<3 



dyj (8) dyk 
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where 

G= ^ Gjkdyj dyk = ^ {ip-,^^dyj®dyk, 

1<«<3 

K = ^ Kjk dyj ®dyk= ^ ~\'d~' d~/ '^^^ ^"^"^ 

2<fc,j<3 2<fcJ<3 

V- r , , 'ST- /dN dN\ , 

L= 2^ Ljkdyj ^dyk = 2^ \-Q—^-Q—/dyj^dyk 

2<fc,i<3 2<k,j<3 

are the first, second and third fundamental forms on dO.. We denote by g^^ the inverse of Qj^. 
Its Taylor expansion, valid in the neighborhood, is given by 

/O \ 

(5^')i<.><3 = /d+ 0(|y|) . (5.6) 

The Lebesgue measure dx transforms into dx = dei^g-jkY/'^dy. This determinant has the Taylor 
expansion 

det{g,kfl^ = l + 0{\y\), (5.7) 
valid in the neighborhood. The magnetic vector potential F = (Fi,F2)-^3) = (—2;2/2, 2:1/2, 0) is 
transformed to a magnetic potential F = (^1,^2,-^3) given as follows, 

F,{y) = j2F,{^-\y))^. (5.8) 
fc=i 

In [22], a particular choice of a gauge transformation is selected so that, in a neighborhood of 
the point xq, the new vector potential F satisfies, 

^1 = 0, F2 = 0(|y|2), F3 = 2/icosz. + y2sinz/ + 0(|2/|2), (5.9) 

Here v = v{xq) is the angle between the magnetic field and the tangent plane of 50 at the point 

Xq. 

Notice that the constants implicit in the O notation in (|5.6p . (|5.7p and (|5.9p can be chosen 
uniform (i.e. independent of the boundary point xq) by compactness and regularity of dVl. 

If u is a function with support in a coordinate neighborhood V ^ we may express the functional 
£q{u,Y) explicitly in the new coordinates as follows, 

ef'iu, F) = / det(5,fe)'/' V g'^dy^ - iKHFj)u x {dy, - inRFk^ 



u 

i<i,fc<3 

2|~|2 I |~|. 
— K \u\ + —\u\ 



2 



dy, (5.10) 



where o $ ^ (with (f) the gauge transformation necessary to pass to the F given 

in dSSD). 

5.3. Decomposition of the energy. Consider a family {xq^i}i of points on the boundary 
50 whose choice will be specified below. For each point xq^i, we may introduce a coordinate 
transformation valid near the point xq,;. Actually, after performing a translation, we may 
reduce to the case corresponding to the coordinate transformation <I> valid near the point xq € 30 
as above. 

We introduce a new partition of unity {xi}i covering the set Oi := suppxij whose support 
will be centered at coordinate neighborhoods of a family of points {xq,;}. Let a = «(«;) <C 1 be 
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a parameter that will be explicitly chosen below. Let, for 6 > 0, 

Os ■■= {(2/1,2/2,2/3) I < 2/1, -5 < 2/2 < 5, -S < ys < 6} . 

We choose {xi}l as smooth non- negative functions such that "^iXii^) = 1 in f^i, Xz = 1 the 
set r2in$;~^(0(i_Q,)5), and suppx; C QiCi^^^ [Os) , and such that there exists a constant C > 
so that 

I 

and 

< C (finite overlap) 

for all X S r^i. 

Using the IMS decomposition formula and the inequality /f^(xK^)^ ~ Xl{x)^)\Xi''P\'^ dx > 0, 
we get the following lower bound of the surface energy, 

£i''iXlh^^J,A) > Y,[£rixihxi^,A)- \Vxl\^\xlh^|;\^ dx^ (5.11) 

Using the bound on Vx/, we bound the error term by 

|2|,, 7,„/.|2 ^ „—2 x—1 



il'lxihil^r dx < . 



We approximate the magnetic potential A by the magnetic potential F. Part of the approxi- 
mation relies on the construction of a suitable gauge transformation. Let 

(l)i{x) = (A(xo,i) - F(xo,/)) • X , 

and 

ui{x) = e*'^'^'^) xi{x)xi{xMx)ilj{x) . 

Then it holds true that, 

S^ixiXihi^, A) = Sl^'im, A - V^i) . (5.12) 
We will prove the following estimate (for large values of k), 

£l''iui,A-V<Pi)>il-5)£l''{ui,F)-C6K''[ (5.13) 

Jn 

The proof of ()5.13p will be postponed to the end of the section. We will bound the energy 
£o{ui,F) from below by expressing it in boundary coordinates. Let z/ = z/; = I'ixQ^i) be the angle 
between the applied magnetic field /3 = (0,0,1) and the tangent plane of dQ at the point xq^i. 
Selecting the gauge giving us the magnetic potential F in (j5.9p . we get with 

ipl = ui o X (a gauge transformation), 

that 

f3°(nz,F)=/" det(9,fc)'/' Yl ^'''idy, - iKHF,)il^i x [dy, - iKHFk)i^i 

J Ox Li ^ • I.^Q 



i<i,'t<3 

dy 



Inserting the estimates (|5.6p and (|5.7p we obtain (again it is assumed that 5 is sufficiently 
small) 

£f'{uuY)>{l-C5)j^ (^|(Vj^_i^i/F(y))^;|'-«Vd' + yl^d') dy 
-C5k^ [ M'^dy. 
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Using the pointwise inequality (with e > arbitrary) 

|(V^ - iKHF)iji\^ > (1 - e)\{Vy - iKHE,^)iPi\^ + (1 - e~^){KHf\{F - E,J^; 
with from (13. 41). we obtain 



£i''iui,F) > il-C£-C6)£l''{i;i,B,^)-Ce-\KHf [ {{F - B,,)^l^l\^ dy 

JOs 



'Os 



We estimate the integral / |(F — F,y^)^pl\'^ dy using ()5.9p . In this way, we find, 

JOs 

e~\KHf I \{F-F,^)i,i\''dy<Ce~\KH)H^ I l^^pdy. 

We conclude that (with the choice e = k6'^) 

El^'iuuF) > {1 - C6^K - C6)£l''{i;i,B,^) - {6 + 6\)k^ [ dy . 

After a scaling, y = {kH)~^/'^ z, we obtain 

£f'{i:u^.,) = ^ [ (|(V.-iE,,(z))^,|2-^|V;zP + -^|V^,|') dz, (5.14) 

where 'ipi{z) = 'iPi[{kH)~^^'^ z^ 

We impose the condition V kH6 ^ 1. Next, by defining A = [k/H — 1]+, we get < k/H < 
1 + A. 

Let b = min [k,/H, 1). It is easy to see that, 

[ ^\iii\^dz>-^= [ b\i^i\^dz-C\{KH) [ \i;i\^dz. 
Jo^, ^ VtiH Jo^, Jos 

We insert this estimate into (|5.14|) , then we apply0 Theorem 13.91 (with i = \J kH6) . In this way, 
we conclude that, 

£i^iui,F)>{l-C5^n-CS)V^E{b,i^{xo^i)){i5^)-C{6 + S'^n + X)K^ [ iV'ipdy. (5.15) 

Jos 

provided that k is large enough. We combine the estimates in ()5.1ip - (|5.15p to get, 

I 
I 

-C{6 + 6^k + X)k^ [ liJil^dy] 
JOs ^ 

We estimate as before using the finite overlap of the supports of the partition of unity, 

V / \ipi\'^dy<C [ \tlj\'^dx<C5. (5.16) 
p JOs Jni 



We need to apply Theorem 13.91 when b £ [Go, 1]. If b G (0, Go), we use Remark 13.51 
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(5.17) 



Next, we note that we have a Riemann sum, 

V|£;(b, z.(xo,0)(45')| - / Eib, v{x)) da{x) < C5 . 
I ^ ^ JDndn 

Hence, we find that 

Sl^ixih^, A) >(1 - C5^K - C6)V^ [_ E{b, i/(x)) da{x) 

JDndn 

- CSV^ - C{S^K + 6 + X)5k^ - Ca~^5~^. 

We choose the parameters 

6 = K-^/\ a = K-^/^^ (5.18) 
so that the remainder terms in (|5.17p are all estimated as 

o^max(K, [k — -ff] + )^ , 

as K — 7- oo. 

Proof of (j5.13p . Let A = A — ^(pi. Using the pointwise inequality 

\zi-Z2\^ > (1 - e)\zi\^ + (1 - e-^)\z2\\ 

which is valid for all complex numbers zi, Z2 and real numbers e £ (0, 1), we obtain the lower 
bound 

£^^iui,A)>{l-e)£l^{ui,F) + e-\KH)^[ \(k - Y)ui\^ dx - e f k^\ui\^ dx. (5.19) 
Observing that 

(A - F)(x) = (A - F)(x) - (A - F){xo,i) , 

we get from Corollarv 14.41 

(KHf [ \ {A-F)ui\'^dx <CX^^^K^ [ \x - xo^il"^ luil"^ dx 
Jn J 



Inserting this into (j5.19p . and making for simplicity the non-optimal choice e = 6, we finish the 
proof of (j5.13p by observing that A < 1 for large values of k. □ 

5.4. The bulk term. We return to the energy decomposition in (j5.4p . The choice we made for 
the parameter 6 allows us to estimate the upper bound in (|5.5p as 0{k^^^) = o{k). The surface 
term in ()5.4p is estimated using (|5.17p . It only remains to estimates the bulk term £Q{x2htp, A) 
appearing in ()5.4p . The estimate of this term was the objective of the first part of this paper 
|12j . Actually, Theorem 6.1 of |12j tells us that, 



£o{X2hi^,A) > E2\D\ [k-H]1 + o{ max(K, [k - H]l)) . 



Inserting this estimate into (j5.4p , then using the lower bound in (|5.17p and the choice of 5 finishes 
the proof of Theorem 15.11 
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6. Upper bound 

The aim of this section is to give an asymptotic upper bound of the ground state energy in 
()1.4p . This will be done through the computation of the energy of relevant test configurations, 
whose construction hints at the expected behavior of the actual minimizers of the energy in (jl.ip . 

The main theorem of this section is stated below. 

Theorem 6.1. For all a{n) = o(l) as k —t- oo and C > 0, there exists < err(K) = o(l) such 
that if 

1 - d{K) < — <C. 

K 

Then, as k — )■ oo, the ground state energy in ()1.4p satisfies, 

Eg,,t{K, H) < / E{b, v{x)) da{x) + E2\9.\[k - H\\ + err ( max(K, [k - H]\)) . 

J do. 

Here b = min [n/H, l), and da{x) is the surface measure on the boundary ofVt. 



Proof of Theorem \6.1[ 
Boundary trial configuration 

Let 5 > he small but fixed. We will choose another parameter r/ > which will be specified as 
a negative power of k below. 

Choose a finite collection of points {xj} C dil such that 

yj k : 5/2 < dist(xj,Xfc) and Vj : mindist(xj, Xfc) < 25. 

Define Uj as 

Uj = {x £ do, : y k ^ j , dist(2;, Xj) < dist(2;, x^)} ■ 

Clearly the Uj^s are disjoint and dQ = [Jj Uj. 

Next, we construct a family of sets that covers a tubular neighborhood of dO. That will 
be done by using the boundary coordinates ?/2, 2/3) introduced in Sec. 15.2.11 {yi = defines 
the corresponding part in dO). Let ^j be the coordinate transformation that straightens a 
neighborhood Vj of the point Xj such that ^j{xj) = 0. We may assume that Uj C Vj for all j. 
Let 

Oj = {x = $-^(yi,y2,y3) : '^J^{0,y2,y3) G Uj and < yi < r?} . 
We now choose a parameter 5': 

^ <^5' <^6 



[5' will be chosen below as a negative power of k). Define Of = ^j{Oj) n {y G : yi = 0}. 
We may cover Of by a square lattice {Kj^i} where each Kj^i is centered at point yj^i and has 
side-length 25' . Let 

Jj = {i : Kj-, C of}, Nj = Card J,- , 
and N = Nj. Clearly, the number N satisfies, 

N X {25' f \dQ\ as (5' ^ . (6.1) 

We combine the coordinate transformation $j by a translation so that the new coordinates 
of the point yj^i are 0. Thus, we let ^j^i be the resulting coordinate transformation valid in 
#J^(M+ X Kj^i) n Vj such that <^j,i{yj,i) = 0. 

We consider only indices i £ Jj for some j. Let Xj^i = $"^^(0). At each point Xj^i, the 
magnetic field (5 = (0,0,1) forms an angle Vj^i = v{xj^i) € [0, 7r/2] with the tangent plane to 
do. As explained earlier in Sec. 15.2. H there exists a real valued function (pj^i such that, if F is 
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the vector field defined in cartesian coordinates by F{xi, X2, X3) = (— 2:2/2, xi/2, 0), and F is the 
vector field defined in y-coordinates by the relation in ()5.8p . then, 

F + V(t)j,i = E^^_^ + Kj^i , 

where 'E^. - is the magnetic potential from (|3.4p , and j is a vector field given in y-coordinates 
by, 

R,,= 0(|yi|2 + |y2|2 + |y3|2) 
\ 0(|yi|2 + |y2|2 + |y3|2) 

Let b = mm{K/H, 1), i = 6W kH and uj^i a minimizer of the functional Q\},Vji,i introduced in 
dsn]). For X = $-/(yi,y2,y3) G Oj, we put, 

0, else 



Notice that by construction the supports of the tpj^i do not overlap. 
We define a test function V'bnd £ H'^{^;C) as follows, 

VxGJ^, V','f(x) = /.(^^^^^^)V'(x), (6.2) 

where /i is a cut-off function satisfying, 

supp/i C [-1,1] , 0</i<linlR, = 1 in [-1/2, 1/2] , 

and 

Energy of the test- configuration 

We will compute the energy of the configuration (^J^^'^, F). Notice that the construction of ijj^^^ 
and the change of variable formulas in Sec. 15.2. II together imply the existence of a constant C > 
(independent of 5) such that, if k is sufficiently large and /3 is an arbitrary real number in (0, 1), 
then we have the upper bound, 

£^^{i,'^f,F)<{l + C{5 + r^ + P))Y, I \{V - inHEu,:)h{yi/r,)uj,,{y^t dy 



' + ] dx , 



where Vs> = (0,oo) x (-5', 6') x (-(5',5')- 
By Theorem 13.61 



(6.3) 



{\uj4z)\^ + \uj4z)\^) dz<Cf, 

where C is a universal constant and £ = 6'^/kH. We use the change of variable formulae in 
Sec. 15.2.11 to express 

iV'jl^ dx + ^ J \ipj\'^dx 

in boundary coordinates, then we apply the change of variable z = y\J kH and use the decay of 
Uj^i. In this way, we get a constant C such that, for all j and i ^ Jj, we have, 

-k" j + y j \ij,,i\Ux<V^ (^-^\u,4z)\' + ^\u,,,\^) dx + -^{6+v){6'f. 
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Here = (0, oo) x {-1,1) x {-1,1). Also, 



n 



Using again Theorem 13.61 

/ |(V - iKHB^^J h{yi/ri) Uj4yV^)\'^ dy 
JVs, 

<(l + /3) / |(V-mi^E,^.Jn,■,(2/^/^)|2dy + /3-V' / [h' {yihi)]^ \uj,i{y^)\^ dy 



< (1 + P)VkH I |(V - iE^^ Juj-il^ dz + I3-^7]-'\kH)-'' I zi\uj,i\' dz 



<{1 + P)^kH / |(V-iE^^Jnj-i|^dz + /3-^r?-^(Ki/)-i(y)2 



Substituting the above upper bounds into (|6.3p . we get, 



if' 



1 



Recall that b = mm{K/H, 1). By our choice of Uj^i, we get that. 



\{V-iB,Juj4' - -\uj^iY + i^\u,,i\ ] dz < (250 [E{h,Vj^^) + o{l)] 



2H 



Also, by (|6.ip . we have 



^{6r < c. 



Inserting this and the bounds on H/k into (|6.4p . we get 

£'''{^P^f,F) < [l + C(5 + r? + /3)]^/^5](25')2[i?(b,^,,) + o(l)] 

+ C[K^r\S'^ + V^) + K{6 + r, + P) + (3-^ri~\~^] . 
We can now, for example, choose 

The sum is a Riemann sum, hence as 5' — )• 0, 

y2i'^S'?E{b,i^j^i) = [ E{b,iy{x))da{x)+o{l). 

Therefore, we get the following upper bound, 

£^^{ip^f, F) < (1 + C6)V^ [ E{b, iy{x)) da{x) + o(k) + C5k . 

Jan 



(6.4) 



(6.5) 



(6.6) 



Since 5 > was arbitrary this is consistent with the boundary part of the upper bound in 
Theorem 16.11 
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Bulk trial configuration 

We keep the choice of the parameters S and i] introduced in the preceding section. Let be 
the test function defined in |12| Eq. (6.15)]. The function V'Jj'^ vanishes in 

{x : dist(2;, dU) < r]} . 

Thus, and ip^^^ have disjoint support. Consequently, by defining, 

/(x) = v:;!^(x)+vj;,f (x), xgo, 

we get. 

Theorem 6.5 in |12) tells us that 

£^^{iP^%,F) < E2\n\ [k - H]l + o(^max (k, [k - H]l)^ . 

Inserting this estimate and that in (|6.6p into (|6.7p . we get, 

'f^°(/,F) < \/^ / E{b,v{x)) da{x) + E2\n\[K - H]\ + o{ m&i^{K,[K - H]\) \ . 

Jdn \ J 

Recalling the ground state energy in (|1.4p . we deduce the following upper bound, 

Eg,st{K, H) < \^ / E{b,iy{x))da{x) + E2\n\[K- H]l + o( mcix{K,[K- H]l)] . 
This finishes the proof of Theorem 16.11 □ 

7. Proof of main theorems 

7.1. Proof of Theorem 11.11 We combine the lower bound of Theorem 15.11 (with D = and 
h = 1) and the upper bound of Theorem 16.11 



7.2. Proof of Theorem 11.21 Let D C 0, he open and smooth. Let fi : M+ — M be a function 
satisfying limK->.oo t-t-ii-^) = 0. Suppose that the magnetic field H satisfies H > k — ^(k)k. 
Using Theorem 3.3 in |12j . we know that, 

Ml^{uj^) = o{l) , asK-^oo, (7.1) 

where 

ujfi = {x£fl : dist(2;, (}il) > (7i(k)/k} , 
and gi{K) is any function satisfying limK^oo5i(^) = 00. 

Lemma 7.1. //(^, A) is a solution of (|1.3p . then, 

£{,l^,A;D)-8{f^,A-D)- [ |V/|Vl'rf^ + v / i^-f)M^dx 

Jd ^ JD 

= -Re[ li^l'^J iy-Vfda + o{K). (7.2) 

holds true as k —)■ 00. Here v is the unit inward normal vector of dO, and f G H^{Q) is any 
function such that, 

V/gL°°(17), supp/C^. 

Proof. Using the estimate in (|7.ip . the result of the lemma follows through an integration by 
parts applied to the term £{f'ip, A; D). See |13| Lemma 6.1] for details. □ 

Lemma 7.2. // {ip, A) is a minimizer of the functional in p.l|) . then as k ^ 00, 

k^H'^ I I curl(A - Y)\^dx = o(max(ft;, [k - H]\)) . 
Jn 
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Proof. Recall the functional £q in (|5.2p . Theorem 15.11 tells us that, 

£:o(V',A)> / E{b,u{x))da{x)+E2\n\[K- H]l + o{m&x{K,[K-H]l)) . 
Jan 

Consequently, the result of the lemma follows by observing that 

£{'iP,A)=£o{^,A) + k'^H^ [ |curl(A-F)|2(ix 

and using the upper bound obtained in Theorem 16.11 □ 

Proof of (jl.Sp . Let ("0) A) be a solution of (|1.3p . Multiplying the first equation in ()1.3p by ■0 
then integrating over D we get, 

/ \{V -iKHA.)^l:\^ dx+ f 4^ u ■ (y - iKH)i; da{x) = [ {1 - dx . (7.3) 

Jd JdD Jd 

Using ()7.ip . we can show that (see \13\ Lemma 6.1] for details), 

ip h' • {V — iKH)ip da{x) = o{k) 



idD 

as K — 7- oo. Consequently, we may rewrite (|7.3p as follows, 

,.2 



£:o(V',A;Z5) = -^ / |V'|^dx + o(K), (7.4) 

where £q is the functional introduced in (|5.2p . 
Let h G C^(M) be a cut-off function such that 

supp/i C [-1,1], /i(x) = 1 in [-1/2, 1/2] , 0</i<linM. 

We define a function / as follows, 

Vx.l., fix) = l-h{^^^^^), (7.5) 

where = \ is the complement in of D, and 

L = max(«;, [k - H]\) . (7.6) 

We will prove the estimate below, 

|V/|Vl''^^+^ / {'^-f?\M''dx+Ke [ \,p\^J u-Vf da = o{ma^{K,[K-H]l)) . (7.7) 
D ^ JD Jan 

Details concerning the derivation of the estimate in (j7.7p will be postponed to the end of this 
proof. 

Using (j7.7p and the decomposition formula of Lemma |7.H we get, 

£o{tl^, A; D) = Soifi^, A; D) + o( max{K, [k - H]l)) . (7.8) 
We apply Theorem 15.11 to bound £Q{f'ip, A; D) from below. In this way, we infer from ()7.8p . 

£oiip,A;D) > \/^ / E{b,iy{x)) da{x) + E2\D\[k - H]l + o{max{n,[K - H]l)) . (7.9) 
JDnan 

Inserting this lower bound into (|7.4p finishes the proof of (jl.Sp . The only point left is the 
justification of the estimate in ()7.7p . 
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Proof of (|7.7p .' Let 51 (k) = L^/^. Recall the estimate in ()7.ip valid in the set w^. In order to 

estimate the integral / (1 - f i.. we use the simple decompoeitio,, 

Jd 

[ {l-f'mUx= I / {l-f')M'dx. 

We estimate the integral over H D as follows: 

/ (1 - /2) VI' dx<o( [ {1- dx) < o{L-'). (7.10) 

It is left to estimate the integral over o;^ n D. Here we use that the measure of the set 
supp(l -/^)n<ni:> C {x G : dist(x,aj^) < c/i(k)/k, dist(x,L»'^) < 2L~^} 

is 



As a consequence, we get by using the bound ||'i/'||L°°(r2) ^ 1 together with our choice of 51 (k) 
that, 

/ (X-f'f\^tdx = 0[L-'l^^-'). 
The estimates obtained for the integrals over H D and uj'^^r\ D together yield that. 



ID 

In a similar fashion, we may show that. 



/ {l-f^f\i^\Ux = 0{K). (7.11) 



f \Vf\M''dx = o{L), Re/ f v ■ V f da = o{L) . 
JD Jdn 

The two aforementioned estimates above, together with that in (j7.1ip and the definition of 
L = max(K, [k — H]'^), yield the estimate in (|7.7p . □ 

Proof of (jl.9p . Suppose that {^,A) is a minimizer of the functional in (jl.ip . Let the function / 
be as in (|7.5p . Since {ip,A) is a solution of (jl.3p . then the formulas in (|7.4p and (|7.8p still hold 
true. 

We define a function g £ H^{Q,) as follows. 



yx^n, g{x) = Vi - /2(x) , 

so that /^(x) + (7'^(x) = 1 in 0. Consequently, we have the decomposition formula, 

f (V^, A) > £:(/V, A) + ^:(5^, A) - / (|V/|2 + I V<7|') IVI' dx . (7.12) 

Is results from Lemma [4. II that as k — t- 00, 

(|V/|2 + |V5|')|^|' dx = o(max(K, [k - H]l)) . (7.13) 

Also, using Theorem 15.11 we obtain the lower bound, 

f(5V,A) > /_ ^(b,z^(2;))(ia(x)+^2|^'|[K-i?]+ + o(max(K, [K-i?]^)). 

We insert this estimate together with that in (|7.13p into (|7.12p . Also, we use the upper bound 
for £{^l^, A) obtained in Theorem 16.11 In this way we infer from (|7.12p . 

£{fi),A)<^^ f_ E{b,v{x))da{x)+E2\D\[n-H]\+o{^&i^{n,[K-H]\)). (7.14) 

J DndO. 
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Using the estimate of Lemma 17.21 we get further, 

Soif^p, A) < [_ E{b, i^{x)) da{x) + E2\D\ [k - H]l + o( max(K, [k - H]l)) . 

JDndn 

We insert this estimate into (|7.8p . This gives us an upper bound of £o{tl^,A;D), which when 
inserted into (j7.4p yields the upper bound, 

-■ ^ / \i^\^dx<V^ [_ E{b,uix))da{x) + E2\D\[k- H]l + o{max{K,[K- H]l)) . 

2 Jd JDndn 

Combining this estimate with (|1.8p finishes the proof of (|1.9p . □ 
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